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1. Introduction 

, Theoretical study of neutron star (NS) cooling began 
with pioneering works of Tsuruta & Cameron (1965) and 
Bahcall & Wolf (1965). It has been argued that one- 
nucleon direct Urea (DU) processes n — *■ peD, pe — )■ ni> 

, are forbidden up to sufficiently high density and the main 
role is played by the two- nucleon modified Urea (MU) pro- 
cesses, like nn — > npeP and np — > ppev. The "standard" 
scenario of NS cooling emerged, where the main process 
responsible for the cooling is the modified Urea process 
MU nn — > npev calculated using free one pion exchange 
(FOPE) between nucleons, see Friman & Maxwell (1979). 
An order of magnitude less contribution to the neutrino 
emissivity is given by the nucleon bremsstrahlung (NB) 
processes. The main process among NB processes respon- 
sible for the cooling is neutron-neutron bremsstrahlung 
nn — > nnvv, neutron-proton bremsstrahlung pn — > pnvv 
is less effective, and the proton-proton one pp ppvv is 
even less effective. This scenario explains only the group 
of slow cooling data. To explain the group of rapidly cool- 
ing data the "standard" scenario was supplemented by 
"exotic" processes either with a pion condensate, with 
a kaon condensate, with hyperons, or involving DU re- 
actions, see Tsuruta (1979), Shapiro & Teukolsky (1983) 
and references therein. All these processes may occur only 
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for densities higher than a critical density, (2 6) np, de- 
pending on the model, where no is the nuclear saturation 
density. 

Pair breaking and formation (PBF) processes per- 
mitted in nucleon superfluids have been suggested. 
Flowers et al. (1976) calculated the emissivity of 15*0 neu- 
tron pair breaking and formation (nPBF) process and 
Voskresensky & Senatorov (1987) considered the general 
case. Neutron and proton pair breaking and formation 
processes (pPBF) were incorporated within a closed di- 
agram technique including medium effects. Numerical 
estimates are valid both for 15*0 and 3P2 superfluids. 
Schaab et al. (1997) have shown that the inclusion of PBF 
processes into the cooling code may allow one to describe 
the " intermediate cooling" group of data (even if one ar- 
tificially suppresses the effects of the medium). Thus the 
"intermediate cooling" scenario arose. 

The PBF processes were then incorporated in the 
cooling codes of other groups, that worked with the "stan- 
dard plus exotics" scenario, see Tsuruta et al. (2002), 
Yakovlev et al. (2004a), Page et al. (2004). Recently 
Page et al. (2004) called the approach that incorporates 
the PFB processes into the "standard" scenario, the 
"minimal cooling" paradigm. Some papers included the 
possibility of internal heating that results in a slowing 
down of the cooling of old pulsars, see Tsuruta (2004) 
and refs. therein. However, because of the simplification 
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of the consideration, calculations within the "standard 
plus exotics" scenario or within the "minimal cooling" 
paradigm did not incorporate effects of the medium. 

The necessity to include in-medium effects in the NS 
cooling problem is a rather obvious issue. It is based 
on requirements of condensed matter physics, of the 
physics of the atomic nucleus and heavy ion collisions, 
see Migdal et al. (1990), Rapp & Wambach (1994), 
Ivanov et al. (2001). The relevance of in-medium ef- 
fects for the NS cooling problem has been shown 
by Voskresensky & Senatorov (1984), (1986), (1987), 
Sonatorov & Voskresensky (1987), Migdal et al. (1990), 
Voskresensky (2001) who calculated emissivity of the MU, 
NB and PFB processes taking into account in-medium 
effects. We call these processes medium-modified Urea 
process (MMU), medium nucleon bremsstrahlung (MNB) 
and medium pair breaking and formation (MPBF) 
processes. The efficiency of the developed "nuclear 
medium cooling" scenario for the description of NS 
cooling was demonstrated within the cooling code by 
Schaab ct al. (1997) and then by Blaschke et al. (2004). 
In the latter paper it was shown that it is possible 
to fit the whole set of cooling data available today. 
Besides the incorporation of in-medium effects into 
the pion propagator and the vertices, it was also ex- 
ploited that the 3P2 neutron gaps are dramatically 
suppressed. The latter assumption was motivated by 
the analysis of the data (see Figs 12, 15, 20 - 23 of 
Blaschke et al. (2004)) and by recent calculations of the 
3P2 neutron gaps by Schwenk & Friman (2004). However 
more recent work of Khodcl ct al. (2004) suggested that 
the 3P2 neutron pairing gap should be dramatically 
enhanced, as a consequence of the strong softening of 
the pion propagator. Thus the results of calculations 
of Schwenk & Friman (2004) and Khodel et al. (2004), 
which both aim to include medium effects in the 
evaluation of the 3P2 neutron gaps, are contradictory. 

Our aim here is to check the consequences of an 
enhanced 3P2 neutron pairing gap within the "nu- 
clear medium cooling" scenario following the work of 
Blaschke et al. (2004). The paper is organized as follows. 
In section 2 we start with a brief recapitulation of the 
Landau-Migdal Fermi liquid approach. Subsection 2.1 cal- 
culates the NN interaction amplitude. In subsection 2.2 
we separate the in-medium pion mode yielding the main 
contribution to the NN interaction for n> uq. Subsection 
2.3 demonstrates the importance of the renormalization of 
the weak interaction in the medium. In subsection 2.4 we 
show internal inconsistencies of the FOPE model, which is 
the basis of the " standard" scenario. Then in section 3 we 
recapitulate the main processes as they arc treated within 
the "nuclear medium cooling" scenario. In section 4 we 
discuss calculations of the pairing gaps and how gaps af- 
fect the emissivities of different processes. Then in section 
5.1 we discuss the cooling model of Blaschke et al. (2004). 
In section 6 we present emissivities of the main processes, 
the heat capacity and thermal conductivity contributions 
in the scenario of Khodel et al. (2004), when the neutron 



processes are assumed to be frozen. In section 7 we show 
our numerical results. Concluding remarks are given in 
section 8. 

2. Medium effects. Nuclear Fermi liquid 

2.1. NN interaction. Hard and soft modes 

At temperatures of interest (T ^ epm Spp) neutrons and 
even protons are only slightly excited above their Fermi 
seas and all the processes occur in the close vicinity of 
the Fermi energies epmSFp- Quasiparticlc approximation 
is fulfilled for nucleons. In such a situation a Fermi liquid 
approach seems to be the most efficient one. Within this 
approach the long-scale phenomena are treated explicitly 
whereas short-scale ones are described by the local quan- 
tities expressed via phenomenological Landau-Migdal pa- 
rameters. We deal with the baryon densities 7 uq Z n Z 
0.5 no. The value no ~ 0.5m^, where m,r — 140 MeV 
is the pion mass, H = c = 1. At such densities related 
to internucleon distances d ~ 1/n^/^ ^ (0.5 1.3) I/tott 
processes important for the description of the NN interac- 
tion correspond to typical energies and momenta u},k ^ 
few rn-Tr. We call them the long-range processes and treat 
them explicitly. These are nucleon particle-hole processes, 
A isobar-nucleon hole processes (since the typical en- 
ergy of the isobar is of the order of the mass difference 
ruA — niiq ~ 2.1 TOtt) and processes related to the exci- 
tation of the pion (the typical excitation energy is of the 
order of the pion mass). Thus we explicitly present loop 
diagrams which depend strongly on the energy and mo- 
mentum for the uj,k ^ few m,r of our interest. 

Using the above argumentation of Fermi liquid theory 
(see Landau (1956), Migdal (1967), Migdal et al. (1990)) 
the retarded NN interaction amplitude is presented as 
follows (see also Voskresensky (2001) for further details) 



X 



+ 




where 




(2) 



The solid line presents the quasiparticlc Green function 
of the nucleon, the double- line is that of the A isobar.^ 
Thus most long-range diagrams of the particle- hole and A- 
nucleon hole types are explicitly presented in (1). Other 
long-range terms come from the pion. The double-wave 
line in (2) corresponds to the exchange of the free pion 

^ Diagrams with open particle lines have sense 
only if the quasiparticle approximation is valid. 
Otherwise one needs to use the closed diagram tech- 
nique developed by Voskresensky & Senatorov (1987), 
Knoll & Voskresensky (1995), (1996). 
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with inclusion of the contributions of the residual S wave 
ttA^'A'' interaction and ww scattering, i.e. the residual inter- 
action irreducible to the nucleon particle-hole and delta- 
nucleon hole. The latter contributions we have taken into 
account in (1). Thus the full particle-hole, delta-nucleon 
hole and pion irreducible block (first block in (2)) is by its 
construction significantly more local than contributions 
given by explicitly presented graphs. If we decided to cal- 
culate the irreducible block in (2) we would need to intro- 
duce a, u), p exchanges and to treat diagrams with multiple 
pion lines. All these terms are less energy-momentum de- 
pendent for u!,k ^ few than those we treat explicitly. 

Up to now we have not perform any approximations, 
except the quasiparticle approximation for the nucleons. 
Now we perform the main approximation. Using the lo- 
cality argument, one assumes that the Landau-Migdal pa- 
rameters, /„„, fnp in the scalar channel and (7„„, gnp in 
the spin channel, which parameterize the first block in (2), 
are approximately constant. Their values are extracted 
from analysis of experimental data. In a more extended 
approach these parameters should certainly be calculated 
as functions of the density, neutron and proton concentra- 
tions, energy and momentum. 

The part of interaction involving A isobar is analo- 
gously constructed 



(3) 



The main part of the NA interaction is due to the 
pion exchange. Although information on the local part 
of the NA interaction is rather scarce, one can con- 
clude (Migdal et al. (1990), Suzuki et al. (1999)) that the 
corresponding Landau-Migdal parameters are essentially 
smaller then those for the NN interaction. Therefore for 
simplicity we neglect the first graph in the r.h.s. of (3). 

Straightforward resummation of (1) in the neu- 
tral channel yields (Voskresensky & Senatorov (1987), 
Migdal et al. (1990)) 




■ a/3 
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= Co {T^p + Z^^a, ■ a^) + flNT^0{<Tr ■ fe)(<72 • fe), (4) 
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Ipp = (1 - 4.gCoA^„)7(.g„„), 



fnn — fpp — f ~^ f •> fnp — / / ; 9nn — 9pp — 9 ~^ 9 1 

and gnp = 9 — 9', the dimensional normalization factor is 

usually taken to be Cq = t^"^ /[mj^ppf^ino)] ~ 300 MeV- 
fm^ ~ 0.77m~^, PFNino) is the nucleon Fermi momentum 



in the atomic nucleus, D^o is the full retarded Green 
function of 7r°, Aaf3 is the corresponding NN~^ loop 
(without spin degeneracy factor 2) 
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Ann ^ -TO*Pi^„(27r2)-i, for UJ < qVFn,q < 2p_F«, PFn = 
rrinVFrn m* is the effective neutron mass in medium, and 
for simplicity we neglect proton hole contributions due to 
the low concentration of protons. The term proportional 
to Co in (4) demonstrates the change of the local block by 
the long-range loop NN correlation factors and the term 
proportional to the pion-nucleon coupling constant /^at 
demonstrates the change of the second (pion) term in (2). 
Resummation of (1) in the charged channel yields 
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+ ZnpO- 1 
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+ fln^npi'^'i- 



fe)(tT2 • k) , 



2/'7(/'), Z^^ = 2g'W), 



T^p = l\9')D^- 



(7) 



(8) 



1 - ^xCoA^p , 



D_ is the full retarded Green function of tt . The 
Landau-Migdal parameters are poorly known for isospin 
asymmetric nuclear matter and for n > no- Taking into 
account the arguments that they can be considered ap- 
proximately as constants, for estimates we may use the 
values extracted from atomic nucleus experiments. One 
can expect that the most uncertain will be the value of the 
scalar constant / due to the; essential role of the medium- 
heavy cr meson in this channel. But this parameter does 
not enter the tensor force channel which is the most im- 
portant in our case. Unfortunately, there are also uncer- 
tainties in the numerical values of the Landau-Migdal pa- 
rameters in other channels and even for atomic nuclei. 
These uncertainties arc, mainly, due to attempts to get 
the best fit of experimental data in each concrete case, 
slightly modifying the parameterization used for the resid- 
ual part of the NN interaction. For example, calcula- 
tions by Migdal (1967) gave / ~ 0.25, /' ~ 1, 5 ~ 
0.5, gi' ~ 1 whereas Saperstein & Tolokonnikov (1998), 
Fayans & Zawischa (1995), Borzov et al (1984), including 
quasiparticle renormalization pre-factors, derived the val- 
ues / ~ 0, /' ~ 0.5 ^ 0.6, g ~ 0.05 ± 0.1, g' ~ 1.1 ± 0.1. 

Typical energies and momenta entering the NN in- 
teraction amplitude of interest are w ~ and k ~ PFn- 
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Then one estimates ^{gnm^ ~ 0, fc ~ ppmn = uq) ~ 
0.35 0.45. For uo — k ~T typical for the weak processes 
with vv one has ^~^{gnm ijJ — k ~T,n — uq) ~ 0.8 0.9. 

2.2. Virtual pion mode 

The pion is strongly affected by the nucleon surroundings. 
It is obvious already from the fact that typical densities are 
n ~ 1 and the pion-nucleon coupling constant is /t^n — 1 
in units = h — c = 1. Obviously perturbation theory 
is not applicable in this case. Instead we continue to use 
the Fermi liquid approach suitable for the description of 
the strong interaction. 

A straightforward resummation of diagrams (1), (2) 
yields the following re-summed Dyson equation for pions 




+ [nil 



(9) 



The ttNA full-dot-vertex includes a phenomenological 
background correction due to the presence of the higher 
resonances, 11^^^ is the residual retarded pion self-energy 
that includes the contribution of all the diagrams which 
are not presented explicitly in (9), such as S wave ttNN 
and TTTT scatterings (included by double- wavy line in (2)). 
For zero temperature the main contribution to fl^g^ is 
given by the Weinberg- Tomozawa term which has a sim- 
ple analytic form. A part of 11^^^ related to tttt fluctu- 
ations is important for the description of the vicinity 
of the pion condensation critical point. It is calculated 
explicitly. Other contributions are rather small and can 
be neglected, as follows from the comparison of the the- 
ory predictions with different atomic nucleus data, see 
Migdal et al. (1990). 

The full ttNN vertex takes into account NN correla- 
tions 





.(10) 



Therefore the nucleon particle- hole part of 11,^0 is cx j{gnn) 
and the nucleon particle-hole part of n^± is cx 7(3')- The 
value of the NN interaction in the pion channel is deter- 
mined by the full pion propagator at small uj and k ^ ppm 
i.e. by the quantity 

(c.*)2(fc) = -(Z?^)-i(c. = 0,fc,/i.). (11) 

Typical momenta of interest are fc ~ pi?„. Indeed the mo- 
menta entering the NN interaction in MU and MMU pro- 
cesses are k = ppn, the momenta governing the MNB 
are k = km (Voskresensky & Senatorov (1986)) where the 
value k = fc„i ~ (0.9 l)pFn corresponds to the mini- 
mum of {ijj*)^{k). The quantity u>* = uj*{km) is the effec- 
tive pion gap, i.e., the effective pion mass in some sense. 
The effective pion gap lo* demonstrates how much the vir- 
tual (particle-hole) mode with pion quantum numbers is 



softened at a given density. The quantity uj*'^ {ppnin)) re- 
places the value (mj -I- p^r^ in the case of the free pion 
propagator used for the calculation of the MU process by 
Friman & Maxwell (1979). It is different for tt^ and for 
TT^ since neutral and charged channels are characterized 
by different diagrams permitted by charge conservation, 
thus also depending on the value of the pion chemical po- 
tential, ^Tj+ 7^ /Itt- 7^ 0, ^7r0 = 0. For T <^ SFn,£Fp, One 
has /i^- = /ig = £i?„ — Spp, as follows from equilibrium 
conditions for the reactions n —^ pi:" and n — > pei>. 

As follows from numerical estimates of different 7 fac- 
tors entering (4) and (7), the main contribution to NN 
interaction for n > no is given by the resummed medium 
one pion exchange (MOPE) diagram 



X 



(12) 



if this channel (T cx {cri ■ k){(T2 ■ k), see (7), (4)) of the 
reaction is not forbidden or suppressed by specific effects 
such as symmetry, small momentum transfer, etc. 

The density dependence of the effective pion gap 
u* that we use in this work is demonstrated in Fig. 
1 (Fig. 1 of Blaschke et al. (2004)). The curve la in 
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Fig. 1. Nucleon - nucleon correlation factor F and square of 
the effective pion gap ui* with pion condensation (branches la, 
2, 3) and without (la, lb). 

Fig.l shows the behavior of the pion gap for n < n^^ , 
where 

is the critical density for the pion condensa- 
tion. In this work for simplicity we do not distinguish 
between different possibilities of 7r°, tt* condensa- 
tions and the so-called alternative-layer-structure, see 
Voskresensky & Senatorov (1984), Migdal et al. (1990) 
and Umeda et al. (1994). Thus we assume that 7r° 
and 7r~ condensations occur at the very same criti- 
cal density n^^ . Although the value n^^ depends on 
rather uncertain different parameters we, following 
Blaschke et al. (2004), further assume n^^ ~ 3 no, cf. 
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discussion by Migdal et al. (1990). The curve lb demon- 
strates the possibihty of a saturation of pion softening 
and the absence of pion condensation for n > n^^ 
(this possibility could be realized, e.g., if Landau-Migdal 
parameters increased with the density). This pion gap 
(from curves la+lb) iw the; value that determines the 
pion Green function for the pion excitations. Curves 2, 
3 demonstrate the possibility of pion condensation for 
n > n^^ . The continuation of branch la for n > n^^' , 
called branch 2, shows the reconstruction of the pion 
dispersion relation in the presence of the condensate 
state. This pion gap is the value that determines the 
pion Green function for the pion excitations above the 
pion condensate vacuum. In the presence of the pion 
condensate (for n > ric-n) the value uj* from curve 2 enters 
the emissivities of all processes with pion excitations in 
initial, intermediate and final reaction states. 

In agreement with the general trend known in con- 
densed matter physics, fluctuations dominate in the vicin- 
ity of the critical point of the phase transition and die out 
far away from it. The; jump from branch la to 3 at n = n^^ 
is due to the first order phase transition to the tt conden- 
sation, see discussion of this point by Dyugaev (1975), 
(1982), Voskresensky & Mishustin (1981), (1982), 
Migdal et al. (1990). The value on branch 3 is 

proportional to the amplitude of the pion condensate 
mean field (the line with the cross in the standard 
notation of the diagram technique). 

The observation that pion condensation appears by the 
first order phase transition needs a comment. The first or- 
der phase transitions in the systems with several charged 
species is associated with the possibility of a mixed phase, 
see Glendenning (1992). The emissivity is increased within 
the mixed phase since efficient DU-like processes due 
to nucleon re-scattering on the new-phase droplets are 
possible. However Voskresensky et al. (2002), (2003) and 
Maruyama et al. (2003), (2005a) demonstrated that, if it 
exists, the mixed phase is probably realized only in a 
narrow density interval due to charge screening effects. 
Thereby to simplify the consideration we further disre- 
gard the possibility of a mixed phase. We also disregard 
the change in the equation of state (EoS) due to pion 
condensation assuming that the phase transition is rather 
weak. 

The density dependence of the correlation factor F is 
presented in Fig. 1. For the correlation factor entering 
the emissivity of the MMU process one can use an ap- 
proximate expression F(n) ~ 1/[1 -I- C(n/no)"'^/^], C ~ 
1.4 ^ 1.6. Note that this value F is an averaged quan- 
tity (F*^ = F^_,rf). Actually the correlation factor F*^ 
entering the emissivity of the MMU process (see eq. (23) 
below) looks more involved and depends on the energy- 
momentum transfer, being different for vertices connected 
to the weak coupling (the correlation factor related to the 
weak coupling vertices F^_s is close to unity) and for ver- 
tices related to the pure strong coupling (Fg is slightly 
less than above introduced factor F), see estimates of the 
energy-momentum dependence of the correlation factors 



at the end of the previous subsection. We see that ver- 
tices are rather strongly suppressed (and this suppression 
increases with the density) but the softening of the pion 
mode is enhanced (w*^ < mj) for n > rid — 0.5 -r- 0.8 no. 
Such a behavior is supported both theoretically and by 
analysis of nuclear experiments, see Migdal et al. (1990). 

Even with full microscopic calculations the functions 
F(n) and w*(n) contain large uncertainties. These un- 
certainties come mainly from simplifications inherent to 
the Landau-Migdal approach to nuclear forces where 
the Landau-Migdal parameters are constants. However it 
seems misleading to disregard medium effects only based 
on existing rather large uncertainties. We believe that fur- 
ther more detailed comparison of theory and experiment 
will allow one to reduce these uncertainties. 

2.3. Renormalization of the weak interaction. 

The full weak coupling vertex that takes into account 
NN correlations is determined by (10) where now the 
wavy line should be replaced by the lepton pair. Thus 
for the vertex of interest, A^i N2lv, we obtain, see 
Voskresensky & Senatorov (1987), Migdal et al. (1990), 

V0 = -^m')k-gAi{9')l<y], (13) 

for the (3 decay and 
G 



Vr 
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pp 



2V2 

G 



[l{fnn)lo - gAl{gnn)l(T] , 



Kpplo — dAlpplc] 



(14) 
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Ipp 



c-v - 2/,ip7(/,m)C'oA 

nn 

(1 - 4gCoAm)7(3rm), 



for processes on the neutral currents N\Ni 

Vr, 



pp 



yN .y-l Q 



(15) 
NIN2VP, 

1.17 • 10-5 GeV-2 is the Fermi 

weak coupling constant, cy = 1 — Asin^Ow, sin^ 9w — 
0.23, gA — 1-26 is the axial-vector coupling constant, and 
Ifj, = u{qi)'yij,{l — 75)^(92) is the lepton current. The pion 
contribution ^ q"^ is small for typical | q |~ T or p^g, and 
for simplicity is omitted. In medium the value of gA (i-e. 
g'^) slightly decreases with the density. 

The 7 factors renormalize the corresponding vacuum 
vertices. These factors are different for different pro- 
cesses. The matrix elements of the neutrino/antineutrino 
scattering processes Nv Nv and of MNB behave 
differently depending on the energy-momentum transfer 
and whether N = n or N = p in the weak coupling 
vertex. Vertices 



V (16) 




N 
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are modified by the correlation factors (5) and (8). 
For N = n these are 'y{gnn,'^,Q) and 7(/nn,<A',9) leading 
to an enhancement of the cross sections for u > qvpn and 
to a suppression for lo < qvpn- 

Renormalization of the proton vertex (vec- 
tor part of Vj^ + Vjjp) is governed by the pro- 
cesses, see Voskresensky & Senatorov (1987), 
Voskresensky et al. (1998), 




being forbidden in vacuum but permitted in medium. For 
the systems with 15*0 proton-proton pairing, oc g\ contri- 
bution to the squared matrix element (see (14)) is com- 
pensated by the corresponding contribution of the dia- 
gram with anomalous Green functions of protons. The 
vector current term is oc Cy in vacuum whereas it is 
oc Kpp in medium (by the first diagram in (17)). Thereby 
the corresponding vertices with protons are enhanced in 
medium compared to small vacuum values (cx Cy ~ 0.006) 
leading to enhancement of the cross sections, up to ~ 
10 -r 10^ times for 1.5 -j- 3no depending on the choice 
of the parameters. It does not contradict the statement 
that correlations are suppressed in the weak interaction 
vertices at n < no- Enhancement with density comes 
from estimation (6) (that directly follows from equation 
(2.30) of Voskresensky & Senatorov (1987)). Also as en- 
hancement factor (up to ^ 10^) comes from the vir- 
tual in- medium photon (7™), second diagram in (17), 
whose propagator contains 1/m^ ^ where rrij 

is the effective spectrum gap, that compensates for 
the small factor from electromagnetic vertices, see 
Voskresensky et al. (1998), Leinson (2000). It is included 
by replacing Vpp. Other processes permitted in 

intermediate states like processes with pp~^ and with 
pions arc suppressed by a low proton density and by 
~ pre-factors, respectively. The first diagram in 
(17) was considered in Voskresensky & Senatorov (1987), 
where the MpPBF process was first introduced, and then 
in Migdal et al. (1990), Schaab et al. (1997), and it was 
shown that MnPBF and MpPBF processes may give con- 
tributions of the same order of magnitude. Finally, with 
electron-electron hole (second diagram (17) and neutron- 
neutron hole (first diagram (17)) correlations included 
we recover this statement and the numerical estimate of 
Voskresensky & Senatorov (1987), Schaab et al. (1997). 



Work of Kolomeitsev & Voskresensky (1999) gives an- 
other example demonstrating that, although the vacuum 
branching ratio of the kaon decays is T{K~ ^ e~ + 
i/e)/T{K~ fj,~ + v^) w 2.5 X 10"'', in medium (due 
to the lambda-proton hole, Ap^^, decays of virtual K~) 
it becomes of the order of unity. Thus we again see that 
depending what reaction channel is considered, in-medium 
effects may significantly enhance the reaction rates or sub- 
stantially suppress them. Ignoring these effects may lead 
to quite misleading results. 

2.4. Inconsistencies of tlie FOPE model 

Since the FOPE model is the basis of the "standard" 
scenario for cooling simulations, which is still in use, we 
would like to demonstrate inconsistencies in the model 
for the description of interactions in dense {n Z no) 
baryon medium (see Voskresensky & Senatorov (1986), 
Voskresensky (2001)). 

The only diagram in the FOPE model which con- 
tributes to the MU and NB is 




Dots symbolize FOPE. This is the first available Born ap- 
proximation diagram, i.e. the second order perturbative 
contribution in /ttat coupling. In order to be theoretically 
consistent one should use perturbation theory up to the 
same second order in /jtat for all the quantities. The pion 
spectrum is then determined by a pion polarization oper- 
ator expanded up to the same order in f^^: 

u'^ c^ml + k^ +U^{u,k,n), (19) 

Il^{w,k,n)= ••( )■■ 

The value Il°{w,k,n) is easily calculated containing no 
uncertain parameters. For w — > 0, A; ~ and for isospin 
symmetric matter 

~ -ao - ipoco, (20) 
^ ^, '^ruNPFk^f^N ^ n o ^. ^IMIn ^ „ 

TT^ TT 

Replacing this value in (19) we obtain a solution ifiow ~ 
{uj*Y{k) with Im w < (for k ~ km) already for n > 
0.3no that would mean appearance of pion condensation 
{{cu*)^{km) < 0). Indeed, the mean field begins to grow as 
~ exp(— Imw • t) ~ exp(ai), a > 0, until repulsive tttt 
interaction does not stop its growth. But it is experimen- 
tally proven that there is no pion condensation in atomic 
nuclei, i.e. even at n = no. 

The puzzle is solved as follows. The FOPE model does 
not work for such densities. One should replace FOPE by 
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the full NN interaction given by (4), (7). The essential 
part of this interaction is due to MOPE with vertices cor- 
rected by NN correlations, see (12). Also the particle-hole, 
NN~^, part of the pion polarization operator is corrected 
by NN correlations. Thus 



(21) 



is suppressed by the factor 7(3', lu = 0, k ^ pp, n ~ no) — 
0.35 ^ 0.45. The final solution of the dispersion relation 
(19), now with full 11 instead of 11°, yields Imo; > for 
n = no whereas the solution with Imcij < 0, which shows 
the beginning of pion condensation, appears only for n > 
n-cTT > no- 



3. Medium effects in neutrino radiation processes 

3.1. Medium effects in two-nucleon processes 

Medium effects essentially modify the contributions of all 
processes. The main contributing diagrams are 



u 




(22) 

It was shown in Voskrcscnsky & Scnatorov (1984), (1986), 
(1987) and Migdal et al. (1990), see also a more recent re- 
view Voskresensky (2001), that the main contribution to 
the MMU process actually comes from the pion channel of 
the reaction nn npeD (the first diagram in (22)), where 
ev are radiated from the intermediate pion exchanging 
nuclcons. Less contribution comes from the NN~^ inter- 
mediate reaction states (second diagram) , and only much 
less contribution for n ^ no comes from the nucleon of the 
leg of the reaction (third diagram, which naturally gener- 
alizes the corresponding MU(FOPE) contribution (18)). 

Moreover, due to the pion softening (medium modifi- 
cation of the pion propagator) the matrix elements of the 
MMU process are further enhanced with the increase of 
the density towards the pion condensation critical point, 
sec Fig. 1. Roughly, the cmissivity of MMU reaction then 
acquires a factor (mainly due to the pion decay channel 



of MMU) 

e.[^M^ 3 .10/3 r6(n) 

where the pre-factor (n/no)^"^*^ arises from the phase 
space volume. 

A different enhancement factor arises for the MNB pro- 
cesses, where radiation from intermediate reaction states 
(see first two diagrams in (22)) is forbidden: 



---t^^NB] . 103 (n/no)^/3 ^"^"^ . 



(24) 



The value lo* entering eqs. (23) and (24) is determined 
by curve la for n < and by curve 2 for n > Uc-k if 
a condensate is present. If a condensate is assumed to be 
absent one should use the continuation lb of the curve la. 



3.2. Medium effects in one-nucleon DU-lil<e processes 

3.2.1. MNPBF processes 

The one-nucleon processes with neutral currents given 
by the second diagram (16) for A'' = {n,p) are for- 
bidden at T > TcN by energy-momentum conservation 
but they can occur at T < T^jv, where T^n is the 
critical temperature for the nucleon-nucleon (NN) pair- 
ing. Then the necessary energy and momentum can be 
taken from breaking and formation of the Cooper pair. 
However they need special techniques to be calculated, see 
Flowers et al. (1976), Voskresensky & Scnatorov (1987). 
Their calculation is easily done in terms of closed dia- 
grams with normal and anomalous Green functions, see 
Voskresensky & Scnatorov (1987). These diagrams con- 
tain only one nucleon loop. It clearly demonstrates 
that these processes are indeed one-nucleon-like processes 
rather than two-nucleon NB-like processes, as one some- 
times interprets them. Due to the one-nucleon origin there 
appears a huge (for the gap A > MeV) pre-factor ~ 10^® 
in their cmissivity, see cq. (37) below. Also the typical tem- 
perature pre-factor ~ following the rough estimate of 
Flowers et al. (1976) is actually misleading. Instead there 
is a A'^(T/A)i/2 pre-factor. 

These processes (MnPBF) n nvv and (MpPBF) 
p — > pv9 play very important roles in the cooling of 
supcrfluid NS, see Voskresensky & Scnatorov (1987), 
Scnatorov & Voskresensky (1987), Migdal et al. (1990), 
Schaab et al. (1997), Blaschke et al. (2004). Due to the 
full vertices in (16) extra r^_, factors appear in the 
MPBF cmissivity. Medium effects are not as important 
for the n — > nvV process, changing the cmissivity by a 
factor of the order of one but they increase the cmissivity 
of the p puD process by two orders of magnitude, see 
eq. (17). 

3.2.2. Pion and kaon condensate processes 

The P wave pion condensate can be of three types: 7r+, 
TT^ , and tt" with different values of the critical densities 
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^cTT = {ncn^ J fi^^^ , TicT^o) , see Migdal (1978). Thus above 
the threshold density for the pion condensation of the 
given type, the neutrino emissivity of the MMU process 
(22) is to be supplemented by the corresponding PU 
processes 




n n 



p n 




The wavy line with the cross is associated with 
the amplitude of the pion condensate mean field. 
Contrary to the FOPE model, the MOPE model of 
Voskresensky & Senatorov (1986) consistently takes into 
account the pion softening effects for n < Uctt and both 
the pion condensation and pion softening effects in pres- 
ence of the condensate for n > Hctt- As we have men- 
tioned, in our numerical calculations we assume for sim- 
plicity that ricTT = n^^ ~ 3 no is the same for neu- 
tral and charged condensates. For vertices in (25) we 
use the values prcscntcid above. Thus, cmissivitics of the 
PU processes are suppressed by the r^r^_g factors. In 
the case of a weak condensate field for non-superfluid 
matter all processes with charged currents yield contri- 
butions to emissivity of the same order of magnitude, 
whereas processes on neutral currents are significantly 
suppressed, see Voskresensky & Senatorov (1984), (1986), 
Leinson (2004). 

For n > n^^ the kaon condensate processes come into 
play. The most popular is the idea of the S wave if ~ con- 
densation (e.g. see Brown et al (1988), Tatsumi (1988)) 
which is allowed at He > i^t-k- (''^k- effective 
kaon mass) due to possibility of the reaction e K~i'. 
Analogous condition for the pion, /Xg > _ (w* _ is the 
effective pion mass) is not fulfilled due to a strong S wave 
TriVA^ repulsion, cf. Migdal (1978), Migdal et al. (1990), 
(again the in- medium effect), otherwise S wave 7r~ con- 
densation would occur at lower densities than K~ con- 
densation. The neutrino emissivity of the K~ condensate 
processes is given by the equation analogous to charged pi- 
ous however with a different NN correlation factor and an 
additional suppression factor due to a small contribution 
of the Cabibbo angle. 

The phase structure of dense NS matter might 
be very rich, including tt*^, tt^ condensates and 
K^, K~ condensates in both S and P waves, 
see Kolomeitsev & Voskresensky (2003); coupling 
of condensates, sec Umeda ct al. (1994); charged 
p-meson condensation, see Voskresensky (1997), 



Kolomeitsev & Voskresensky (2004); fermion conden- 
sation yielding an efficient DU-like process in the vicinity 
of the pion condensation point (with the emissivity 
~ lO^^^ Tg^ mil oc 1/T), see Voskresensky et al. (2000); 
hyperonization, see Takatsuka & Tamagaki (2004); quark 
matter with different phases, such as 2SC, CFL, CSL, 
plus their interaction with meson condensates, see 
Rajagopal & Wilczek (2000), Blaschke et al. (2001), 
Grigorian et al. (2004) and refs therein; and dif- 
ferent mixed phases. In the present work, as in 
Blaschke et al. (2004), we suppress all these possibil- 
ities of extra efficient cooling channels except a PU 
process on the pion condensate. Other choices are 
eff'ectively simulated by our PU choice. 



(25) 3.2.3. Other resonance processes 



There are many other reaction channels allowed in 
the medium. Any Fermi liquid permits propagation of 
zero sound excitations of different symmetry related 
to the pion and the quanta of a more local interaction 
determined via Landau-Migdal parameters fa.fj and 
ga,i3- These excitations present at T ^ may also 
participate in the neutrino reactions. The most essential 
contribution comes from the neutral current processes, 
see Voskresensky & Senatorov (1986), given by the first 
two diagrams of the series 



-I- 



-I- 



.(26) 



Here the dotted line is the zero sound quantum of the 
appropriate symmetry. These are the resonance processes 
(second one of the DU-type) analogous to those processes 
on the condensates with the only diflFerencc that the rates 
of reactions with zero sounds are proportional to the ther- 
mal occupations of the corresponding spectrum branches 
whereas the rates of the condensate processes are pro- 
portional to the modulus squared of the condensate mean 
field. The contribution of the resonance reactions is rather 
small due to the small phase space volume [q ^ T) associ- 
ated with zero sounds. Note the analogy of the processes 
(26) with the corresponding phonon processes in the crust. 

3.2.4. DU processes 

The proper DU processes in matter, as n — > pe~Ve and 




(27) 



p p n 

should also be treated with the full vertices. They are 
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forbidden up to the density n)r when the triangle in- 
equality pFn < PFp + PFe begins to be fulfilled. For tra- 
ditional EoS like that given by the variational theory by 
Akmal et al. (1998), DU processes are permitted only for 
n > 5 no- Due to the full vertices in (27) extra r^_j fac- 
tors appear in the DU emissivity. 

4. Gaps 

In spite of the many calculations that have been per- 
formed, the values of nucleon gaps in dense NS matter 
are poorly known. This is the consequence of the expo- 
nential dependence of the gaps on the density dependent 
potential of the in-medium NN interaction. This potential 
is not sufficiently well known. Gaps that we have adopted 
in the framework of the "nuclear medium cooling" sce- 
nario, see Blaschke et al. (2004), are presented in Fig. 2. 
Thick dashed lines show proton gaps used in the work 
of Yakovlev et al. (2004a) performed in the framework of 
the "standard plus exotics" scenario. In their model, pro- 
ton gaps are artificially enhanced (that is, not supported 
by any microscopic calculations) to get a better fit of the 
data. We use their "Ip" model. Neutron 3P2 gaps pre- 
sented in Fig. 2 (thick dash -dotted lines) are the same as 
those of the "3nt" model of Yakovlev et al. (2004a). We 
will call this model I. Thin lines show 15*0 proton and 3P2 
neutron gaps from Takatsuka & Tamagaki (2004), for the 
model AV18 by Wiringa et al. (1995) (we call it model 11). 
We take the same l^o neutron gap in both models 1 and II 
(thick solid line), as calculated by Ainsworth et al. (1989) 
and used by Schaab et al. (1997) within the cooling code. 
Blaschke et al. (2004) used both models I and 11 within 
the "nuclear medium cooling" scenario. Since the ISq neu- 
tron pairing gap exists only within the crust, dying for 
baryon densities n > 0.6 uq, its effect on the cooling is 
rather minor. The effect on the cooling arising from the 
proton ISq pairing and from the neutron 3P2 pairing, 
with gaps reaching up to rather high densities, is pro- 
nounced. The NS cooling essentially depends on the val- 
ues of the gaps and on their density dependence. Findings 
of Schulze et al. (1996) and Lombardo & Schulze (2000), 
who incorporated in-medium effects, motivated us to 
check the possibility of a suppression of ISq neutron and 
proton gaps. For that we introduced pre-factors for l^o 
neutron and proton gaps which we varied in the range 
0.2 1, see Figs. 18 and 19 of Blaschke et al. (2004). 

4.1. Possibilities eitlier of strong suppression or strong 
enhancement of 3P2 g3p 

Recently Schwenk & Friman (2004) have argued for a 
strong suppression of the 3P2 neutron gaps, down to val- 
ues ^ 10 keV, as a consequence of the medium-induced 
spin-orbit interaction. They included important medium 
effects, such as the modification of the effective interac- 
tion of particles at the Fermi surface owing to polarization 
contributions, with particular attention to spin-dependent 
forces. In addition to the standard spin-spin, tensor and 




Fig. 2. Neutron and proton pairing gaps according to model I 
(thick solid, dashed and dotted lines) and according to model 
II (thin lines), see text. The ISo neutron gap is the same in 
both models, taken from Ainsworth et al. (1989). 



spin-orbit forces, spin non- conserving effective interac- 
tions were induced by screening in the particle-hole chan- 
nels. Furthermore a novel long-wavelength tensor force 
was generated. The polarization contributions were com- 
puted to second order in the low-momentum interaction 
Mowk- These findings motivated Blaschke et al. (2004) to 
suppress values of 3P2 gaps shown in Fig. 2 by an ex- 
tra factor /(3P2, n) = 0.1. Further possible suppression of 
the 3P2 gap is almost not refiected in the behavior of the 
cooling curves. 

Contrary to expectations of 

Schwenk & Friman (2004), the more recent work of 
Khodel et al. (2004) argued that the 3P2 neutron pairing 
gap should be dramatically enhanced, as a consequence of 
the strong softening of the pion propagator. According to 
their estimation, the 3P2 neutron pairing gap is as large 
as 1 10 MeV in a broad region of densities, see Fig. 1 
of their work. 

In order to apply these results to a broad density in- 
terval both models may need further improvement. The 
model of Schwenk & Friman (2004) was developed to de- 
scribe not too high densities. It does not incorporate 
higher order nucleon-nucleon hole loops and the A iso- 
bar contributions and thus it may only partially include 
the pion softening effect at densities uq. The model of 
Khodel et al. (2004) uses a simplified analytic expression 
for the effective pion gap (ti;*)^(fcm), valid near the pion 
condensation critical point, if the latter occurred by a sec- 
ond order phase transition. The latter assumption means 
that (ci;*)^(fcm) is assumed to be zero in the critical point 
of the phase transition. Outside the vicinity of the critical 
point the parameterization of the effective pion gap that 
was used can be considered only as a rough interpolation. 
Actually the phase transition is of first order and evalu- 
ations of quantum fluctuations done by Dyugaev (1982) 
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show that the value of the jump of the effective pion gap 
in the critical point is not as small. Moreover repulsive 
correlation contributions to the NN amplitude have been 
disregarded. In the pairing channel under consideration, 
already outside the narrow vicinity of the pion condensa- 
tion critical point, the repulsion originating from the NN 
correlation effects may exceed the attraction originating 
from the pion softening. Notice that, if the pairing gap 
enhancement occurred only in a rather narrow vicinity of 
the pion condensation critical point, it would not affect 
the results of Blaschke et al. (2004). In the latter work 
two possibilities were considered (see Fig. 1): i) a satu- 
ration of the pion softening with increase of the baryon 
density resulting in the absence of pion condensation and 
ii) a stronger pion softening stimulating the occurrence of 
the pion condensation for n > ric — 3 no- In both cases 
the effective pion gap was assumed never to approach zero 
and undergoing a not too small jump at the critical point 
from a finite positive value ((w*)^ — 0.3 m^) to a finite 
negative value ((w*)^ — —0.1 ttv^). The reason for such a 
strong jump is the strong coupling. If it were so, the strong 
softening assumed by Khodel et al. (2004) would not be 
realized. However, due to uncertainties in the knowledge 
of forces acting in strong interacting nuclear matter and 
a poor description of the vicinity of the phase transition 
point we cannot exclude that the alternative possibility 
of a tiny jump of the pion gap exists. Therefore we will 
check how these alternative hypotheses may work within 
our "nuclear medium cooling" scenario. Avoiding further 
discussion of the theoretical background of the models, we 
investigate the possibility of a significantly enhanced 3P2 
neutron pairing gap and of a partially suppressed proton 
15*0 gap, as suggested by Khodel et al. (2004). To proceed 
in the framework of our "nuclear medium cooling" sce- 
nario we introduce the enhancement factor of the original 
3P2 neutron pairing gap /(3P2, n), and a suppression fac- 
tor of the proton l^o gap f{lSo,p)- We do not change the 
neutron l^o gap since its effect on the cooling is minor. 

4.2. Suppression of the emissivity by superfluid gaps 

Generally speaking, the suppression factors of superfluid 
processes are given by complicated integrals. As was 
demonstrated by Sedrakian (2005) with the example of 
the DU process, these integrals are not reduced to the 
i?-factors, sec Yakovlev et al. (2004a). However, for tem- 
peratures significantly below the critical temperature the 
problem is simplified. With an exponential accuracy the 
suppression of the specific heat is governed by the factor 
^nn for neutrons and for protons: 

^u^exp[^AuiT)/T], for T kT^; i = n,p, (28) 

and = 1 for T > T^, Td is the corresponding critical 

temperature. Wc do not need higher accuracy to demon- 
strate our result. Therefore we will use these simplified 
factors. 

For the emissivity of the DU process the sup- 
pression factor is given by min{^„„, ^pj,}, see 



Lattimer et al. (1991). The same suppression factor 
appears for other one nucleon processes, as PU and KU 
and the resonance processes on zero sounds, see the sec- 
ond diagram of (26). Suppression factors for two nucleon 
processes follow from this fact and from the diagrammatic 
representation of different processes within the closed 
diagram technique by Voskresensky & Scnatorov (1987) 
and Knoll k Voskresensky (1995, 1996). These are: 
^nn • min{^„„, ^pp} for the neutron branch of the MU 
process (and for the same branch of the medium mod- 
ified Urea process, MMU); ^pp • min{^„„, ^pp} for the 
corresponding proton branch of the process; for 
the neutron branch of the (medium modified) nucleon 
bremsstrahlung (MnB) and ^pp for the corresponding 
proton branch of the bremsstrahlung (MpB). Thus, for 
Ann ^ App both neutron and proton branches of the 
MMU process are frozen for T <^ Ten due to the factors 
and ^pp^nni respectively. Zero sound and other 
phonon processes shown by the first diagram (26) are not 
suppressed by the factors. However their contribution 
to the emissivity is very small due to the smallness of the 
available phase space volume. 

5. Cooling model of Blaschke et al. (2004) 

5.1. EoS and structure of NS interior, crust, surface 
5.1.1. NS interior 

We will exploit the EoS of Akmal et al. (1998) (specif- 
ically the Argonne V'18 + 6v + UIX* model), which 
is based on the most recent models for the NN in- 
teraction with the inclusion of a parameterized three- 
body force and relativistic boost corrections. Actually we 
adopt a simple analytic parameterization of this model 
given by Hcisclberg & Hjorth- Jensen (1999) (HHJ). It 
uses the compressional part with the compressibility K ~ 
240 MeV, and a symmetry energy fitted to the data 
around the nuclear saturation density, and smoothly in- 
corporates causality at high densities. The density depen- 
dence of the symmetry energy is very important since 
it determines the value of the threshold density for the 
DU process. The HHJ EoS fits the symmetry energy 
to the original Argonne V18 + 6v + UIX* model yield- 
ing r;.;^u ~ 5.19 „o {M^^ ~ 1.839 Mq). The original 
Argonne EoS allows for neutral pion condensation (for 
n > 2no) that only slightly affects the energy density. 
One may disregard this small change. This EoS does not 
allow for charged pion condensation. The HHJ param- 
eterization of the EoS does not include tt condensation 
effects. Wc further assume that pion condensation (neu- 
tral and charged) occurs for n > 3no, see discussion in 
Blaschke et al. (2004). We assume a minor effect of the 
pion condensation on the EoS and disregard it. We also 
disregard changes of the isotopic composition due to the 
charged pion condensation. The latter effect would be 
small only if the charged pion condensate field were rather 
weak. Thus we assume that the HHJ parameterization of 
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the EoS includes both mentioned effects or that they are 
negUgible. 

5.1.2. NS crust 

The density n ~ 0.54-0.7 no is the boundary of the NS in- 
terior and the inner crust. In the latter there occurs the so- 
called "pasta phase" discussed by Ravenhall et al. (1983), 
see also Maruyama et al. (2005b). Then there is the outer 
crust and the envelope. Note that our code generates a 
temperature profile that inhomogeneous during the first 
10^ lO'^ y. The influence of the crust on the cooling and 
heat transport is minor due to its rather low mass con- 
tent. Therefore the temperature also changes slightly in 
the crust up to the envelope. 

5.1.3. Envelope 

The resulting cooling curves depend on the Tin — 7s re- 
lation between internal and surface temperatures in the 
envelope. Fig. 3 shows uncertainties existing in this re- 
lation. A calculation is presented for the canonical NS: 




Iog„(TJK]) 

Fig. 3. The relation between the inner crust temperature 
and the surface temperature for different models. Dash-dotted 
curves indicate boundaries of the uncertainty band. Notations 
of lines are determined in the legend. For more details see 
Blaschke et al. (2004) and Yakovlev et al. (2004b). 

M = 1.4Mq, i? = 10 km with the crust model HZ90 of 
Yakovlev et al. (2004b). Below we will show that a min- 
imal discrepancy of the increased gap scenario with the 
data is obtained with what we called "our fit" model 
in Blaschke et al. (2004). Using other choices like the 
"Tsuruta law" (r^™-^ = (lOTin)^/^, where and Tin are 
measured in K) only increases the discrepancy. To com- 
pare results with "our fit" model we use the upper bound- 
ary curve, 77 = 4 • 10^^^ and the lower boundary curve 
77 = 4 • 10^* . In Fig. 3 we also draw lines rj — 1 ■ 10^^^ 
and T] = 1 ■ 10^^^ as they are indicated in the corre- 
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sponding Fig. 2 of Yakovlev et al. (20G4b). The selection 
of 77 = 4 • 10~^ and 77 = 4 • 10~^® as the boundaries of 
the uncertainty-band seems to be a too strong restriction, 
see Yakovlev et al. (2004b). The limit of the most massive 
helium layer is achieved for rj ^ 10^^°. On the other hand 
the helium layer begins to affect the thermal structure only 
for ?7 > 10~^^. Thus one could exploit 10~^^ < 77 < 10"^" 
as a Tin — Tg band. We will use a broader band, as is shown 
in Fig. 3. By this we simulate the effect of maximum un- 
certainties in the knowledge of the Tin ^ relation. 

5.2. Main cooling regulators 

We compute the NS thermal evolution adopting our fully 
general relativistic evolutionary code. This code was orig- 
inally constructed for the description of hybrid stars by 
Blaschke et al. (2001). The main cooling regulators are 
the thermal conductivity, the heat capacity and the emis- 
sivity. In order to better compare our results with results 
of other groups we try to be as close as possible to their 
inputs for the quantities that we did not calculate our- 
selves. Then we add the changes, improving the EoS and 
including the medium effects. 

5.2.1. Thermal conductivity 

We take the electron-electron contribution to the ther- 
mal conductivity and the electron-proton contribution for 
normal protons from Gnedin & Yakovlev (1995). The to- 
tal contribution related to electrons is then given by 

1/Ke = l/Hee + l/K^p. (29) 

For T > Tcp (normal "n" matter), we have hc^^ — f^ep- 

For 

T < Tcp (superfluid "s" matter) we use the expression 

'*ep — f^ep/£,pp > Kgp, (30) 

that gives a crossover from the non-superfluid case to the 
superfluid case. The vanishing of k^^ for T <C Tcp is a 
consequence of the scattering of superfluid protons on the 
electron impurities, see Blaschke et al. (2001). Following 
(29) we get k" < k|. 

For the neutron contribution, 

Kji = 1/ Knn "i" 1/ ^np i (^1) 

we use the result of Baiko et al. (2001) that includes cor- 
rections due to the superfluidity. Although some medium 
effects are incorporated in this work, the nucleon-hole 
corrections of correlation terms and the modification of 
the tensor force are not included. This should modify 
the result. However, since we did not calculate k„ our- 
selves, we can only roughly estimate the modification: 
not too close to the critical point of the pion condensa- 
tion the squared matrix element of the NN interaction 

l^^lmod ~ PFn^'i/i^*)'^^ ^'^^ (12) and values shown 
in Fig.l, is of the order of the corresponding quantity 
|Af|^ac PF,ri/[™^ '^Ppn] estimated with the free one 
pion exchange, whereas l-Mlj^j^^,^ may significantly increase 
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for n ~ n™. However Blaschke et al. (2004) checked that 
both increasing and decreasing of the thermal conductiv- 
ity does not change the picture as a whole. 

The proton term is calculated similarly to the neutron 

one, 

Kp = 1/ Kpp -\- 1/ K^ip ' (*^^) 

The total thermal conductivity is the sum of the partial 
contributions 

For the values of the gaps used in Blaschke et al. (2004) 
the other contributions to this sum arc smaller than 
those presented explicitly (Kg, Kn and Kp). Finally 
Blaschke et al. (2004) concluded that in their scenario 
transport is relevant only up to the first 10^ y. 

5.2.2. Heat capacity 

The heat capacity contains nucleon, electron, photon, 
phonon and other contributions. The main in-medium 
modification of the miclcon heat capacity is due to the 
density dependence of the effective nucleon mass. We use 
the same expressions as Schaab et al. (1997). The main 
regulators arc the miclcon and the electron contribu- 
tions. For the nucleons {i = n,p), the specific heat is 
(Maxwell (1979)) 

a ~ 102°(m*/miv) {rii/noy/^^ii Tg erg cm-^K"! , (34) 
for the electrons it is 



5.2.3. Emissivity 

Wc adopt the same set of partial cmissivitics as in the 
work of Schaab ct al. (1997). The phouou contribution to 
the emissivity of the 3P2 supcrfluid phase, as well as the 
zero sound contribution, is negligible. The main emissivity 
regulators arc the MMU, see the rough estimation (23), 
MnPBF and MpPBF, and MNB processes. For n > n™ ~ 
3 no the PU process becomes efficient and for n > nj?^ ~ 
5.19 no the DU process is the dominant. 

Only the qualitative behavior of the interaction shown 
in Fig. 1 is motivated by the microscopic analysis whereas 
the actual numerical values of the correlation parameter 
and the pion gap are rather uncertain. Thus we vary the 
values r(n) and w*^(n) in accordance with the discussion 
of Fig. 1. By that we check the relevance of alternative 
possibilities: a) no pion condensation and a saturation of 
the pion softening with increasing density, curves la-l-lb, 
and b) the presence of pion condensation, curves la+2-|-3. 
We also add the contribution of the DU for n > n]?^. 

All cmissivitics are corrected by correlation effects. 
The PU process contains an extra P^ factor compared to 
the DU process (the emissivity of the latter is oc P^_^). 
Another suppression of PU emissivity comes from the fact 
that it is proportional to the squared pion condensate 
mean field Near the critical point ^ 0.1, increas- 
ing with density up to jy'p ~ /^/^i where f„ ~ 93 MeV 
is the pion decay constant. Finally, the PU emissivity is 
suppressed by about 1-2 orders of magnitude compared 
to the DU one. We adopt the same gap dependence for 
the PU process as for the corresponding DU process. In 
supcrfluid matter all emissivity terms are suppressed by 
the corresponding factors. 



Ce ~ 6 X 10^^ (ne/no)^/^ Tg erg cm"^ K" 



(35) 



Near the phase transition point the heat capacity ac- 
quires a pion fluctuation contribution. For the first order 
pion condensation phase transition this additional contri- 
bution contains no singularity, like for the second order 
phase transition, see Voskresensky & Mishustin (1982) 
and Migdal et al. (1990). Finally, the nucleon contribution 
to the heat capacity may increase up to several times in 
the vicinity of the pion condensation point. The effect of 
this correction on global cooling properties is unimportant 
and simplifying we neglect it. 

The symmetry of the 3P2 supcrfluid phase allows for 
the contribution of Goldstone bosons (phonons): 



i3 K' 



(36) 



for T < Tc„(3P2), n > ncn(3P2)- We include this term 
in our study although its effect on the cooling is minor. A 
similar contribution comes from other resonance processes 
permitted on zero sounds, see the first diagram in (26). In 
order not to introduce an extra parameter dependence we 
will simulate the effect of all phonon and zero sound terms 
with the term (36). 



6. Main cooling regulators in the scenario of 
Khodel et al. (2004) 

When neutron processes are frozen the most efficient pro- 
cess is the MpPBF process, p pvv, for T < Tcp. 
Taking into account the medium effects in the weak cou- 
pling vertex, see eq. (17), we use the same expression 
for the emissivity of this process as has been used by 
Voskresensky (2001) and Blaschke et al. (2004): 



£^ [MpPBF] r 

rj, -,1/2 



10 



29 





PFp 






rUN 


.PFn{no)_ 




MeV_ 



n 7 



■'PP 



2 

■PP 



erg 



cm'^ sec 



T <T, 



cp' 



(37) 



This process contributes only below the critical tem- 
perature for proton pairing. Inclusion of medium effects 
greatly enhances the vertex of this process compared to 
the vacuum vertex, see the above discussion of this fact in 
subsection 2.3. A factor r^_5 ^ 10^ arises, since the pro- 
cess may occur through nn~^ and ee~^ correlation states, 
with subsequent production of vP from the nn~^uP and 
ee^^vv channels rather than from a strongly suppressed 
channel pp~^i'v, see Voskresensky & Senatorov (1987), 



H. Grigorian and D.N. Voskresensky: Medium effects in cooling of neutron stars and the 3P2 neutron gap 



13 



Senatorov & Voskresensky (1987), Migdal et al. (1990), 
Voskresensky et al. (1998), Leinson (2000) and 
Voskresensky (2001). Relativistic corrections incor- 
porated in the description of the pp~^vD vertex in 
Yakovlev et al. (2004a) also produce an enhancement 
being however significantly less than that arises from 
medium effects in nn~^vi> and ee~^vv channels. Thus we 
see no reason not to include these medium effects and 
we note only a moderate dependence of the result on the 
uncertainties in the strong interaction. 

We also present here an explicit expression for 
the emissivity of the proton branch of the nucleon 
bremsstrahlung including medium effects, MpB, pp 
ppvv. In the case of suppressed neutron 3P2 gaps 
this process contributes much less than the nucleon 
bremsstrahlung processes involving neutrons. However, 
when neutron reactions are frozen, the pp — > ppvD pro- 
cess becomes the dominating process for Ten > T > Tcp. 
The emissivity of the pp — > ppvi? reaction takes the form, 
see Voskresensky & Senatorov (1986) for more details, 



Hi] (IL) 



pp pp{oj*Y{km) 

5/3 



erg 



cm-^ sec 



(38) 



Tg — T/10^ K, m* is the effective proton mass. Here we 
take Ts ~ T{n) ~ 1/[1 -I- C{n/noY'% C ~ 1.4 1.6. 
This factor takes into account NN correlations in strong 
interaction vertices, Yp = Up/n is the proton to nucleon 
ratio. As above, for simplicity we assumed that the value 
k = km (at which the square of the effective pion gap 
(w*)^(A:) is minimimum) is rather close to the value of 
the neutron Fermi momentum ppn, as it follows from the 
microscopic analysis of Migdal et al. (1990). To simplify 
the consideration we take the same value of the effective 
pion gap for the given process as that for the MMU process 
(although in general it is not so, and thus the result (38) 
is model dependent), cf. Blaschke et al. (2004), 



64 \PFpJ PFn 



(39) 



We have checked that for T < Tcp for the pairing gaps 
under consideration the MpB reaction contributes signifi- 
cantly less than the MpPBF process. It could be not the 
case only in a narrow vicinity of the pion condensation 
critical point, if pion condensation occurred with only a 
tiny jump of the effective pion gap in the critical point. 
However, even in this case there are many effects which 
could mask this abnormal enhancement. 

For n > n^^ the process pir^ — > pi^v on the neutral 
pion condensate is still permitted. However its contribu- 
tion is strongly suppressed, see Leinson (2004). 

In the case of frozen neutron degrees of freedom the 
specific heat is governed by protons and electrons, see eq. 
(34) for i = p and (35). Here, we again suppress the con- 
tribution to the specific heat by the narrow vicinity of the 



pion condensation critical point due to the fact that in our 
scenario (see Fig. 1) the modulus of the squared effective 
pion gap (w*)^ is always larger than ~ (0.1 0.3) mj. 
With such an effective pion gap the pion contribution to 
the specific heat is not as strong and can be disregarded 
to simplify the consideration. For the second order phase 
transition (and for a first order phase transition but with a 
tiny jump of |(w*)^| in the critical point), pion fluctuations 
would contribute more strongly to the specific heat yield- 
ing a term cx T^^ /lu* for \T - T^^\/T^^ < T^^ , see 
Voskresensky & Mishustin (1982), Migdal et al. (1990). 

When neutron processes are frozen the values k„ and 
Kp are suppressed and the thermal conductivity is reduced 
to the electron and proton contributions, see eqs. (29) and 
(32). 

7. Numerical results 

Now we give the results of our calculations of cooling 
curves. First we present Fig. 21 of Blaschke et al. (2004), 
now Fig. 4. Cooling curves shown in this figure were cal- 





:J0822 







































1 .250 — 
1,320 
1.400 



Fig. 4. Fig. 21 of Blaschke et al. (2004). Gaps are from Fig. 2 
for model II. The original 3P2 neutron pairing gap is addition- 
ally suppressed by the factor f{3P2,n) = 0.1. The pion gap is 
determined by curves la-|-2-|-3 of Fig. 1. The Ts — Tin relation 
is given by "our fit" curve of Fig. 3. Here and in all subsequent 
figures the value Ts is the red-shifted temperature. NS masses 
are indicated in the legend; see Blaschke et al. (2004). 

culated using "our fit" model of the crust, shown with 
the solid curve in Fig. 3. Here and in the correspond- 
ing figures below the surface temperature is assumed to 
be red-shifted, as observed at infinity from the radia- 
tion spectrum. Gaps are given by model II of Fig. 2. 
However, the 3P2 gap is additionally suppressed by a 
factor f{3P2,n) = 0.1, as motivated by calculations of 
Schwenk & Friman (2004). The calculation includes pion 
condensation for n > n™. The figure shows a good fit 
to the data. If we used calculations disregarding the pos- 
sibility of pion condensation (see curves la-|-lb of Fig. 
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1) we would also get an appropriate fit to the data, of. 
Fig. 20 of Blaschke et al. (2004). If we took the original 
3P2 gap of model II, we would not be able to describe 
the data. The cooling would be too fast, see Fig. 23 of 
Blaschke et al. (2004). 

We will now check the possibility of ultra-high 3P2 neu- 
tron pairing gaps, as proposed by Khodel et al. (2004). In 
Figs. 5 and 6 we demonstrate the sensitivity of the results 
presented in Fig. 4 to the enhancement of the neutron 3P2 
gap and to a suppression of the ISq proton gap, following 
the suggestion of Khodel et al. (2004). Again we use the 
calculation including pion condensation for n > n^^ . 




Fig. 5. Cooling curves according to the nuclear medium cool- 
ing scenario, see Fig. 4. Gaps are from Fig. 2 for model II but 
the 3P2 neutron pairing gap is additionally enhanced by a fac- 
tor /(3P2,n) — 50 and the ISo proton gap is suppressed by 
f{lSo,p) = 0.1. The pion gap is determined by curves la-|-2-|-3 
of Fig. 1. The Tg — Tin relation is given by "our fit" curve of 
Fig. 3. 

We start with "our crust" model and model II for 
the gaps, using however the additional enhancement fac- 
tor f(SP2,n) — 50 for the neutron 3P2 gap. Introducing 
factors /(lS'o,p) = 0.1 and f{lSo,p) = 0.5 we test the 
sensitivity of the results to the variation of the 15*0 proton 
gap. We do not change the value of the 15*0 neutron gap 
since its variation almost does not influence on the cool- 
ing curves for NSs with masses M > 1 Mq, that we will 
consider. 

Comparison of Figs. 4-6 shows that in all cases NSs 
with masses M > 1.8 Mq cool in similar ways in spite 
of the fact that 3P2 neutron and 15*0 proton gaps vary 
over wide limits. This is because 3P2 neutron and ISo 
proton gaps disappear at the high densities occuring in 
the central regions of these very massive NSs, see Fig. 
2. Thus these objects cool similarly to non-superfluid ob- 
jects. Extremely rapid cooling of stars with M > 1.84 AIq 
is due to the DU process, which is very efficient in normal 
matter. The latter process appears in the central region 




Fig. 6. Same as Fig. 5, but for the original ISo proton gap 
suppressed by f{lSo,p) ~ 0.5. 



of NSs with M > 1.839 Mq. Therefore we notice that 
the cooling curves are very sensitive to the density depen- 
dence of the gaps. The difference in the cooling of NSs 
with M < 1.8 Mq in the cases presented in Figs. 5 and 6 
is the consequence of different values of proton gaps used 
in these two calculations. This difference is mainly due 
to the MpPBF processes. The larger the proton gap, the 
more rapid is the cooling. 

We checked that for stars with M ^ 1.6 Mq for T < 
Ten for the 3P2 neutron pairing, a complete freezing of 
the neutron degrees of freedom occurs for /(3P2, n) > 20. 
Then contributions to the emissivity and to the specific 
heat involving neutrons are fully suppressed. For heavier 
stars (M > 1.6 Mq) a weak dependence on the value of 
the factor /(3P2, n) still remains even for /(3P2, n) > 100 
but the corresponding cooling curves lie too low to allow 
for an appropriate fit of the data. This difference between 
cooling of stars with M < 1.6 Mq and M > (1.6-^-1.7) Mq 
is due to the mentioned density dependence of the neutron 
3P2 gap. The latter value smoothly decreases with the 
increase of the density reaching zero for n ^ 4.5 no (the 
density 4.5 no is achieved in the center of a NS of mass 
M — 1.7 Mq). At densities slightly below 4.5 uq the gap 
is small. Therefore for stars with M > (1.6-^ 1.7) Mq the 
scaling of the gap by a factor /(3P2,n) changes the size 
of the region where the gaps may affect the cooling. For 
stars with M ^ 1.6 Mq gaps have finite values even at the 
center of the star. Thus there exists a critical value of the 
factor /(3P2,n), such that for higher values of /(3P2,n) 
the cooling curves are unaffected by its change. 

Figs. 5 and 6 demonstrate that we did not succeed in 
reaching an appropriate overall agreement with the data; 
the cooling was too rapid. The cooling of old pulsars is 
not explained in all cases. Although the heating mecha- 
nism used by Tsuruta (2004) may partially help in this 
respect, the discrepancy between the curves and the data 



H. Grigorian and D.N. Voskresensky: Medium effects in cooling of neutron stars and the 3P2 neutron gap 



15 



points seems to be too high, especially in Fig. 6. We see 
that in our regime of frozen neutron processes a better fit 
is achieved in Fig. 5, i.e., for a stronger suppressed proton 
gap (for f(lSo,p) — 0.1). The discrepancy is even more 
severe, since to justify the idea of Khodel et al. (2004) we 
should use a softer pion propagator. Only a strong soften- 
ing of the pion mode might be consistent with a significant 
increase of the neutron 3P2 gap- On the other hand such 
an additional softening would result in still more rapid 
cooling. The work of Voskresensky et al. (2000) discussed 
the possibility of a novel very efficient process with the 
emissivity e^, oc T^, that would occur due to non- Fermi 
liquid behavior of the Fermi sea very near the pion con- 
densation critical point with the assumption of strong pion 
softening. If we included this very efficient process, the dis- 
agreement with the data could be strongly enhanced. An 
enhancement of the specific heat due to pion fluctuations 
within the vicinity of the pion condensation point cannot 
compensate for the acceleration of the cooling owing to the 
enhancement of the emissivity. Khodel et al. (2004) used 
the value ric = 2 uq for the critical density of the pion 
condensation. In case of the Urbana-Argonne equation of 
state that we exploit here (we use the HHJ fit of this 
equation of state that removes the causality problem, see 
Blaschke et al. (2004) for details) the density n — 2 uq 
is achieved in the central region of a NS with the mass 
M ~ 0.8 Mq. This means that aU NSs with M > 0.8 M© 
would cool extremely fast and would not be seen in soft 
X rays. 

Actually, we checked the whole interval of variation of 
/(3P2,n) and f{lSo,p) in the range 1 100 and 0.1-^-0.5 
respectively. We verified that the variation of f{3P2,n) 
and f{lSo,p) in the whole mentioned range done within 
our parameterization of the effective pion gap does not 
improve the picture (curves la-|-2-|-3 of Fig. 1). Using 
the pion gap from branches la-flb does not achieve a 
better fit. In all cases we obtain too rapid cooling. To 
demonstrate this, in Fig. 7 we show the cooling of a 
1.4 Mq star for different values of the f{iP2,n) factor. 
The factor f{lSo,p) is taken to be 0.1. We see that for 
f{3P2,n) < 15 20 the curves rise with the increase of 
/(3P2,n). For f{3P2,n) > 20 the curves do not depend 
on/(3P2,n). 

To check how the results are sensitive to uncertainties 
in our knowledge of the value (39) that determines the 
strength of the in-medium effect on the emissivity of the 
MpB process we multiplied (38) by a pre-factor /(MpB) 
that we varied in the range /(MpB) = 0.2 -^ 5. In agree- 
ment with the above discussion, for /(MpB) < 1, for tem- 
peratures logTs[K] > 5.9 the cooling curves are shifted up- 
wards. For /(MpB) > 1, for temperatures logrs[K] > 5.9 
the cooling curves are shifted downwards. However, in- 
dependently of the value /(MpB) for logrs[K] < 5.9 the 
curves are not changed. Thus it does not allow us to di- 
minish the discrepancy with the data. 

Now we will check the efficiency of another choice of 
the gaps, as motivated by model I, thick lines in Fig. 2. 
Compared to model II model I uses an artificially en- 
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Fig. 7. Cooling curves of the neutron star with mass 1.4 Mq 
according to the nuclear medium cooling scenario, see Fig. 4. 
Gaps are from Fig. 2 for model II but the 3P2 neutron pair- 
ing gap is additionaUy enhanced by different factors /(3P2,n) 
(shown in Figure) and the ISo proton gap is suppressed by 
/(15'o,p) ~ 0.1. The pion gap is determined by curves la-|-2-|-3 
of Fig. 1. The Ts — Tin relation is given by "our fit" curve of 
Fig. 3. 

hanced proton gap. Thus, one can expect that model I 
is less realistic than model II. Also note a different den- 
sity dependence of the proton gaps in model I (it cuts 
off at densities n > 3.2 no) and in model II (it cuts off 
at higher densities, n > 4.2 uq). Since uncertainties in 
the existing calculations of the gaps are very high and 
the parameterization of model I has been used by one of 
the groups working on the problem of cooling of NSs, see 
Yakovlev et al. (2004a), we will consider the consequences 
of this possibility as well. 

Fig. 8 demonstrates our previous fit of the data within 
model I, but for the original 3P2 neutron gap suppressed 
by /(3P2,n) = 0.1, see Blaschke et al. (2004). Again the 
pion gap is determined by curves la-|-2-|-3 of Fig. 1. 

If we took the original 3P2 gap of model I, we would 
not succeed in describing the data. The cooling would be 
too fast, see Fig. 22 of Blaschke et al. (2004). Therefore 
we check the possibility of ultra-high 3P2 neutron pairing 
gaps, as proposed by Khodel et al. (2004). 

As in Fig. 5, Fig 9 uses /(3P2,n) = 50 and f{lSo,p) = 
0.1 but now for the gap in model I, and, as in Fig. 6, Fig. 
10 uses /(IS'ojP) = 0.5 for the gap in model I. Figs. 9 and 

10 show that within the variation of the gaps of model I 
the discrepancy to the data is still stronger compared to 
that for the above calculation based on the use of model 
II. The difference between curves shown in Figs. 9 and 10 
is less pronounced than for those curves in Figs. 5 and 6. 
Indeed, as we have mentioned, the density dependence of 
the proton gap is different in models I and II. In model 

11 the proton gap reaches higher densities (~ 4.2 uq) than 
in model I (~ 3.2 no). Thus in the case shown by Figs. 
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Fig. 8. Fig. 15 of Blaschke et al. (2004). Gaps from model I. 
The original 3P2 neutron gap is suppressed by /(3P2, n) — 0.1. 
The pion gap is determined by curves la+2+3 of Fig. 1. The 
Ts — Tin relation is given by "our fit" curve of Fig. 3. For more 
details see Blaschke et al. (2004). 



9 and 10 a non-superfluid core begins to contribute at 
lower values of the star mass. Therefore in both figures 
the corresponding coohng curves are almost the same for 
M > I.GMq. Using the pion gap from branches la+lb, 
i.e. disregarding the possibility of pion condensation, does 
not allow a better fit. 
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Fig. 9. Cooling curves according to the nuclear medium cool- 
ing scenario, see Fig. 8. Gaps are from Fig. 2 for model I but 
the 3P2 neutron pairing gap is additionally enhanced by a fac- 
tor f{3P2,n) — 50 and the l^o proton gap is suppressed by 
f{lSo,p) = 0.1. The pion gap is determined by curves la-|-2-|-3 
of Fig. 1. The Tg — Tin relation is given by "our fit" curve of 
Fig. 3. 



Fig. 10. Same as Fig. 9, but for the original ISo proton gap 
suppressed by f{lSo,p) — 0.5. 



The dependence of the results on the different choices 
of the Ts — Tin relation is demonstrated by Figs. 11 and 
12 for gaps based on a modification of model II. For 
this demonstration we first took the upper boundary 
curve 77 = 4 • 10~^^ and then the lower boundary curve 
7] = 4.0 • 10^* in Fig. 3. We show that these choices how- 
ever do not allow us to improve the fit. Comparing Figs. 
11 (yy = 4 • 10^^^) and 5 ("our fit") based on the same 
modification of model II we see that with the "our fit" 
crust model (Fig. 5) the deviation from the data points is 
less pronounced. We have checked that based on model I 
one arrives at the same conclusion. An increase of the fac- 
tor f{lSo,p) to 0.5 reduces the fit. In Fig. 12 we use the 
lower boundary curve 77 = 4.0 • 10~* of Fig. 3. We further 
demonstrate that the selection of a different choice of the 
Ts — Tin relation within the band shown in Fig. 3 does not 
diminish the discrepancy. This discrepancy increases com- 
pared to that demonstrated by "our fit" model. Indeed, 
the cooling evolution for times t ^ 10^ yr (Tg lO^K) is 
governed by neutrino processes. Thus the higher Tin, the 
larger Tg is. The slowest cooling is then obtained, if one 
uses the lower boundary curve ry = 4.0 • 10~* of Fig. 3. 
The evolution of NSs for times t ^ 10^ yr begins to be 
controlled by the photon processes. In the photon epoch 
{t ^ 10^ yr) the smaller the T, value, the less efficient the 
radiation is. Thus for t ^ 10^ yr the slowest cooling is ob- 
tained if one uses the upper boundary curve rj = 4.0-10"^^ 
of Fig. 3. The "our crust" curve simulates the transition 
from one limiting curve to the other demonstrating the 
slowest cooling in the temperature interval shown in the 
figures. 

In all cases the data are not explained by the assump- 
tion of an enhanced neutron SP2 gap (for f(3P2,n) > 1) 
and a partially suppressed ISq proton gap (for f{lSo,p) = 
0.1^0.5;. 



H. Grigorian and D.N. Voskresensky: Medium effects in cooling of neutron stars and the 3P2 neutron gap 



17 




Fig. 11. Cooling curves according to the nuclear medium cool- 
ing scenario. Gaps are from Fig. 2 for model II but the 
3P2 neutron pairing gap is additionally enhanced by a fac- 
tor f{3P2,n) — 50 and the ISo proton gap is suppressed by 
f{lSo,p) = 0.1. The pion gap is determined by curves la-|-2-|-3 
of Fig. 1. The Ta — Tin relation is given by the crust model for 
r; = 4.0 ■ 10"^^ 




Fig. 12. Same as in Fig. 11, but for the crust model rj = 4.0 • 
10"**. 

8. Concluding remarks 

Our aim was to consider large 3P2 gaps within the "nu- 
clear medium cooling" scenario of Blaschke et al. (2004) 
that well described the cooling data with the oppo- 
site assumption of suppressed 3P2 gaps. Therefore we 
did not incorporate extra assumptions of internal heat- 
ing for old pulsars, see Tsuruta (2004), and quark 
cores in massive NSs, see Blaschke et al. (2001) and 
Grigorian et al. (2004) and refs therein. 

The main problem with the given scenario is that 
at the frozen neutron contribution to the specific heat 
and to the emissivity, the region of surface temperatures 



Ts > lO^K is determined by proton processes. The most 
efficient among them is the MpPBF process for T < Tcp 
and MpB for T > Tcp. For the proton gaps that we deal 
with, the MpPBF process proves to be too efficient, yield- 
ing too rapid cooling. Thus at least several slow cooling 
data points (at least in data for old pulsars) are not ex- 
plained. Some works ignore the medium-induced enhance- 
ment of the MpPBF emissivity that results in a 10-100 
times suppression of the rate. We omitted this possibility 
as physically unrealistic. The origin of this enhancement 
is associated with opening up of new reaction channels in 
the medium that are forbidden in vacuum. 

Thus, we have shown that the "nuclear medium cool- 
ing" scenario of Blaschke et al. (2004) fails to appropri- 
ately fit the neutron star cooling data with the assump- 
tion of a strong enhancement of the 3P2 neutron gaps (we 
checked the range /(3P2, n) — 1-^ 100) and for moderately 
suppressed 15o proton gaps (for f{lSo,p) = 0.1 0.5). 
On the other hand the same scenario allowed us to ap- 
propriately fit the whole set of data with the assump- 
tion of a significantly suppressed 3P2 neutron gap (for 
f{3P2,n) ~ 0.1). We observed essential dependence of the 
results not only on the values of the gaps but also on their 
density dependence. We used the density dependence of 
the gaps according to models I and II. The latter model is 
supported by microscopic calculations. We excluded an at- 
tempt to artificially fit the density dependence of the gaps 
trying to match cooling curves with the data. Although 
such an attempt could improve the fit, it is not physical 
and we did not pursue it. However we encourage further 
attempts of microscopic calculations of the gaps, which 
would take into account the most important medium ef- 
fects. With correctly treated gaps one could perform new 
simulations of NS cooling. 
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